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The quasiparticle effective mass is a key quantity in the physics of electron gases, describing the 
renormalization of the electron mass due to electron-electron interactions. Two-dimensional electron 
gases are of fundamental importance in semiconductor physics, and there have been numerous 
experimental and theoretical attempts to determine the quasiparticle effective mass in these systems. 
In this work we report quantum Monte Carlo results for the quasiparticle effective mass of a two- 
dimensional homogeneous electron gas. Our calculations differ from previous quantum Monte Carlo 
work in that much smaller statistical error bars have been achieved, allowing for an improved 
treatment of finite-size effects. In some cases we have also been able to use larger system sizes than 
previous calculations. 

PACS numbers: 73.20.-r, 71.10.Ay, 02.70.Ss 



I. INTRODUCTION 

Two-dimensional (2D) electron gases are ubiquitous 
in modern semiconductor devices. Surprisingly, how- 
ever, there remain significant gaps in our knowledge of 
the properties exhibited by these fascinating systems. 
In recent years experimentalists have realized increas- 
ingly high-quality (low-disorder) 2D homogeneous elec- 
tron gases (HEGs) at low densities in quantum-well 
structures 1 ^ and field-effect transistors^ As the density 
is lowered, correlation effects play an ever more impor- 
tant role, and it is hoped that the resulting exotic be- 
havior could be exploited in a new generation of elec- 
tronic and spintronic devices. At very low densities the 
Coulomb repulsion between the electrons dominates, and 
the HEG forms a Wigner crystal^— At higher densities 
the kinetic energy dominates and the electrons form a 
Fermi fluid, in which most properties are qualitatively 
(and at high densities, quantitatively) similar to those 
of a free-electron gas. Fermi liquid theory^ is the phe- 
nomenological framework within which the normal be- 
havior of Fermi fluids is understood. 

According to Fermi liquid theory, low-lying excitation 
energies in HEGs are free-electron-like, and the effects 
of interactions are encapsulated in (i) a renormalization 
of the electron mass (the quasiparticle effective mass) 
and (ii) a set of parameters describing the interaction 
of pairs of excited quasiparticles. In this work we use 
quantum Monte Carlo (QMC) methods to determine the 
quasiparticle effective mass by calculating single-particle 
excitation energies as differences in the total energy when 
electrons are either added to or removed from the ground 
state. 

In our calculations we have used the variational 
Monte Carlo (VMC) and diffusion Monte Carlo (DMC) 
methods^ In VMC we take the expectation value of 



the many-electron Hamiltonian with respect to a Slater- 
Jastrow-backflow trial wave function? 1 ^— which is opti- 
mized by minimizing first the variance of the energ y 14 : 15 
then the energy expectation valued with respect to free 
parameters in the wave function. In DMC 17 we simulate 
a population of "walkers" whose dynamics are governed 
by the Schrodinger equation in imaginary time in order 
to project out the ground-state component of an initial 
wave function. The fixed-node approximation^ is used 
to impose fermionic antisymmetry. All our QMC calcu- 
lations were performed using the CASINO code^ 

In Ref. [20| we presented a DMC calculation of the 2D 
HEG single-particle energy band, enabling us to predict 
the quasiparticle effective mass. In the present work we 
have had access to the Jaguar machine at Oak Ridge 
Leadership Computing Facility, enabling us to achieve 
higher accuracy in our DMC calculations, and leading to 
a refinement of our earlier work. 

The rest of this paper is structured as follows. In Sec. 
HI! we give an overview of the relevant aspects of Fermi 
liquid theory. In Sec. IIIII we describe our computational 
approach. Our results are presented in Sec. IIV1 Finally, 
we draw our conclusions in Sec.[V] We use Hartree atomic 
units, in which the Dirac constant, the electronic charge 
and mass, and 47r times the permittivity of free space are 
unity (h = \e\ = m e = 4ireo = 1), throughout. 



II. LANDAU ENERGY FUNCTIONAL 

A. Parameterization of excitation energies 

According to Fermi liquid theory the total energy of a 
HEG in a particular excited state is given by the Landau 
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energy functional 

E = Eo+^e^SN^a 

k,(T 

+ \ ]T U,*WWit,*Mv,*> (l) 

(k,<7)^(k',(7') 

where JA/k.tr is the change to the ground-state quasipar- 
ticle occupation number for wavevector k and spin a, and 
Eq is the ground-state energy. The energy band £ a (k) is 
the energy of an isolated quasiparticle. Near the Fermi 
surface, the energy band may be assumed to be linear 
and hence we may write 

£ a (k)=£ F + ^(k-k F ), (2) 

TO* 

where £f is the Fermi energy, hp is the Fermi wavevector, 
and to* is the quasiparticle effective mass. The Landau 
interaction function f a%a i (k, k') describes energy contri- 
butions arising from pairs of quasiparticles, and will not 
be considered further in this paper. 

The goal of this work is to obtain accurate values for 
the 2D HEG quasiparticle effective mass to* in the ther- 
modynamic limit at different densities and for different 
spin polarizations, giving us the most important contri- 
bution to the Landau energy functional. 

B. Spin-polarization effects 

Both theoretical wor k 20 i 21 and experimental studies^ 
have shown that the quasiparticle effective mass has a sig- 
nificant dependence on the spin polarization of the HEG. 
We have calculated the effective mass for both param- 
agnetic and ferromagnetic (fully spin-polarized) HEGs. 
Fully spin-polarized HEGs are experimentally relevant 
because they may be created by applying an in-plane 
magnetic field to a 2D electron system. Differences in 
the quasiparticle effective masses of ferromagnetic and 
paramagnetic HEGs result in differences in the trans- 
port properties, which could be exploited in electronic 
or spintronic applications, e.g., in devices that use the 
spin-Coulomb-drag effect^ 

C. Finite-size errors 

The 2D HEGs encountered in real devices are suffi- 
ciently large that they can be regarded as being of essen- 
tially infinite extent. In QMC simulations we can only 
study small numbers of electrons, however. For a HEG 
in a finite simulation cell subject to periodic boundary 
conditions, momentum quantization limits the available 
wavevectors {k} to a discrete lattice. Furthermore, long- 
range Coulomb and correlation effects cannot be treated 
exactly in a finite cell ) 23 i 24 giving rise to finite-size errors 
in the energy band and hence effective mass. 



Fermi liquid theory is only valid for excitations near 
the Fermi surface: in this region the quasiparticle lifetime 
becomes large and hence the quasiparticle momentum oc- 
cupancies are good quantum numbers £ The energy band 
is defined by the Landau energy functional at all k, but 
does not correspond to the quasiparticle band except in 
the vicinity of the Fermi surface. In the infinite-system 
limit, the exact energy band is smooth in general and, if 
the quasiparticle effective mass is well-defined, the band 
must be at least differentiable at the Fermi surface. 



III. QMC CALCULATIONS 

A. Choice of simulation cell 

In all our calculations the simulation cell was square 
and the simulation-cell Bloch vecto r 25 i 26 was k s = 0. 
The number of electrons in the ground state was cho- 
sen to give a closed-shell configuration in each case. For 
ferromagnetic HEGs, our calculations were performed 
with N = 29, 57, and 101 electrons in the ground state. 
For paramagnetic HEGs our calculations were performed 
with ./V = 26, 50, 74, and 114 electrons in the ground 
state. 

The simulation cell was identical for all excitations of 
a given HEG; hence the electron density increased when 
electrons were added and decreased when electrons were 
removed from the ground-state configuration. This pro- 
cedure results in zero finite-size error for a free-electron 
gas. 

B. Trial wave functions 

We use real, single-determinant trial wave functions for 
the closed-shell ground states, which is a computationally 
efficient approach that facilitates the optimization of the 
wave function. In our QMC calculations we used Slatcr- 
Jastrow-backflow trial wave functions. The Jastrow fac- 
tors consisted of polynomial and plane-wave expansions 
in the interelectron distances^ while the backflow func- 
tions consisted of polynomial expansions in the interelec- 
tron distances The polynomial expansions were cut off 
smoothly at the radius of the largest circle that could 
be inscribed in the simulation cell. The Jastrow factor 
and backflow function contained a total of 35 and 17 free 
parameters, respectively, for paramagnetic HEGs, and 
27 and 8 free parameters, respectively, for ferromagnetic 
HEGs. Extrapolation of the VMC energy with different 
trial wave functions to zero energy variance^ suggests 
that our DMC calculations retrieved more than 99% of 
the correlation energy. 

For each density, system size, and spin polarization the 
wave function was optimized in the ground state and the 
resulting Jastrow factor and backflow function were used 
in all the excited states, with the exception of a couple 
of test cases, as discussed in Sec. IIV A 21 
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C. DMC time steps, etc. 

The DMC time steps used in our calculations were 
0.04, 0.2, and 0.4 a.u. at r s = 1, 5, and 10, respectively, 
for paramagnetic HEGs, and 0.01, 0.2, and 0.4 a.u. at 
r s = 1, 5, and 10, respectively, for ferromagnetic HEGs. 
It was verified that halving the time step had a negligible 
effect on the energy band: leading-order time-step errors 
cancel out of the total-energy differences involved. The 
target population exceeded 1200 configurations in each 
case, ensuring that population-control bias is negligible. 

The number of equilibration steps discarded from the 
start of each DMC calculation was sufficiently large that 
the root-mean-square distance diffused by each electron 
in the equilibration period exceeded the linear size of the 
simulation cell. 



IV. RESULTS 
A. Energy band 

1. DMC results for the energy band 

The DMC energy band of the paramagnetic 2D HEG, 
obtained with a Slater- Jastrow-backflow trial wave func- 
tion, is shown in Fig. [TJ Analogous results for a fully 
ferromagnetic HEG are shown in Fig. [3J 

The energy bands calculated in this work at any given 
system size are in agreement with those calculated in Rcf. 
|20| . but the statistical error bars in the present work are 
very much smaller due to the considerably larger com- 
putational resource available. Furthermore, the random 
noise in the trial wave function due to optimization by 
energy minimization with a finite sampling of configu- 
ration space is greatly reduced because of the enormous 
numbers of configurations that were used in the opti- 
mizations. As a result it is now possible to discern a 
systematic trend in the energy band with system size, 
with the bandwidth tending to increase with system size 
N . This in turn leads to a reduction in the predicted 
quasiparticle effective mass in the thermodynamic limit, 
as discussed in Sec. IIVBI 



0.5 



a 

03 

5-0.5 



-1.5 
-0.18 



= -0.20 
§-0.22 



-0.24 



-0.120 



3 



g -0.125?- 



-0.130 



(a) r =1 






Free electron 




Hartree-Fock 


• 


N = 2(, 


■ 


N = 5Q 




N = 74 


A 


N= 114 


-4 


N= 162 





JV= 26 (reopt) 




N = 14 (reopt) 




FIG. 1: (Color online) Energy bands £ (k) for paramagnetic 
2D HEGs of density parameter (a) r s = 1, (b) r B = 5, and 
(c) r s — 10 at different system sizes N. For the curves la- 
beled "reopt," the wave function was optimized separately in 
the ground state and excited states. The free-electron and 
Hartree-Fock bands have been offset so that they coincide 
with the DMC bands at kp. The inset to panel (a) shows the 
energy band around k = in greater detail. 



2. Effect of reoptimizing the wave function in excited states 

The excitation of a single electron or pair of electrons 
has no effect on the optimal Jastrow factor or backflow 
function in the thermodynamic limit; hence the fact that 
the Jastrow factor and backflow function can be reop- 
timized in an excited state in a finite cell is simply a 
manifestation of finite-size error. Reoptimizing the wave 
function when an electron was subtracted from k = in 
a 26-electron HEG at r s = 1 lowered the DMC energy 
by 0.000854(4) a.u., reducing the DMC bandwidth and 
hence increasing the finite-size error, as can be seen for 



N = 26 in the inset of Fig. [lja). In our calculations we 
therefore optimized the trial wave function in the ground 
state and then continued to use the same Jastrow factor 
and backflow function in our excited-state calculations. 
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FIG. 2: (Color online) As Fig.[TJ but for ferromagnetic HEGs. 



B. Quasiparticle effective mass 

1. Quartic fits to the energy bands 

At each density, spin polarization, and system size a 
quartic function £(k) = ao+a2k 2 +oi4k 4 was fitted to the 
DMC energy-band values. The effective mass [defined in 
Eq. ([2])] was then calculated as 



ki 



(d£/dk) kf 



(3) 



tive mass as a function of the range Afc about the Fermi 
wave vector over which we perform the fit. The figures 
also show the effective mass when energy-band data from 
within ±10% of ftp are excluded from the fit. It is clear 
that the effective mass becomes pathological when Afc 
becomes small (i.e., only excitations in the vicinity of fcp 
are considered), for the reasons discussed briefly in Sec. 
IIV B 21 and at length in Ref. HU As the width of the re- 
gion over which the fit is performed becomes larger, the 
effective mass estimates settle down to well-defined val- 
ues that behave in a systematic fashion with system size. 
There is no evidence of any need to exclude data from 
around kp however: effectively, the fitting process aver- 
ages out the pathological behavior in the vicinity of the 
Fermi surface. 



2. Extrapolation to the thermodynamic limit 

The quasiparticle effective masses are plotted against 
system size in Fig. [5] A systematic trend in the effective 
mass as a function of system size can be seen. Hence 
we are able to extrapolate the effective mass to the ther- 
modynamic limit, significantly reducing finite-size errors. 
The effective-mass values that we report in this work are 
exp ected to be more accurate than those reported in Ref. 

m 

The finite-size error does not show the iV -1 / 4 behavior 
predicted by Holzmann et al££ for excitations near the 
Fermi surface, presumably because we have fitted to the 
entire band. (Any TV -1 / 4 behavior in the band near the 
Fermi surface is masked by the pathological behavior that 
is seen in Hartree-Fock theory and hence QMC<2£) To in- 
vestigate the behavior of the effective mass as a function 
of system size, we have fitted the function 



m*{N) = m*(oo) + WV~ 



(4) 



We have investigated the dependence of the estimate 
of the effective mass on the range of energy-band data 
used to perform the fit. Figures [3] and @] show the effec- 



where m*(oo), b, and 7 are fitting parameters, to the raw 
data shown in Fig.|6]and we have performed repeated fits 
to the data with Monte Carlo sampling of the error bars. 
We find that 7 = 1.8(4), 1.4(4), and 1.4(3) at r s = 1, 5, 
and 10, respectively This indicates that the optimal ex- 
ponent 7 is between 1 and 2. This conclusion is reinforced 
by the results shown in Table U where we examine the 
X 2 values of the fits and the extrapolated effective masses 
when different exponents are used in Eq. The extrap- 
olation shown in Fig. [5] assumes an exponent of 7 = 3/2. 



3. Quasiparticle effective mass as a function of density 

Results for the 2D HEG quasiparticle effective mass 
obtained by different authors, including the present work, 
are shown in Figs. [6] and [7] for paramagnetic and fully 
ferromagnetic HEGs, respectively. The present results 
are also given in Table |Hl For ferromagnetic HEGs 
our revised effective masses are in reasonable agreement 
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FIG. 3: (Color online) Quasiparticle effective mass m* against 
range of wavevectors included in the quartic fit of the en- 
ergy band. Specifically, only wavevectors in the interval 
[Isf — Afc, kF + Afc] are used in the fit. For the dashed lines, 
wavevectors within 10% of kF are excluded from the fit. 
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FIG. 4: (Color online) As Fig.O but for ferromagnetic HEGs. 
For N = 57 electrons, excluding wavevectors within 10% of 
kF eliminates all the wavevectors above fcf at which energy- 
band data are available. 



TABLE I: Extrapolated quasiparticle effective mass m*(co) 
and x 2 value for fits to the effective-mass data as a function of 
system size for paramagnetic HEGs with different exponents 
7 in the finite-size fitting formula [Eq. Q]. 



7 


m" 


(oo) (a.u 


) 




X 2 (a.u.) 




r s = 1 


r B — 5 r a 


= 10 


r s = 1 


r s — 5 r 


, = 10 


1/4 


0.88 


0.3 - 


0.7 


21.5 


9.6 


12.4 


1/2 


0.924 


0.74 


0.3 


16.8 


5.8 


7.7 


1 


0.947 


0.97 


0.85 


9.1 


1.11 


1.52 


3/2 


0.955 


1.04 


1.03 


5.2 


0.28 


0.32 


2 


0.959 


1.08 


1.13 


5.3 


2.03 


3.1 



with our previously published results^ For paramag- 
netic HEGs, however, the finite-size errors are relatively 
large at low density, and hence finite-size extrapolation 
reduces the effective masses at r s = 5 and 10 by a sig- 
nificant amount. For paramagnetic HEGs we now find 
that the effective mass remains close to 1 (i.e., electron- 
electron interactions result in almost no renormalization 
of the electron mass) at all the r s values we have consid- 
ered. 

Once again we emphasize that there is no significant 
disagreement between the data reported in the present 
article and Ref. HO- The revision of the effective mass 
simply results from the fact that the random noise in our 
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FIG. 5: (Color online) Quasiparticle effective mass m* against 
N~ 3/2 , where N is the system size, for 2D HEGs. 



TABLE II: Quasiparticle effective masses for paramagnetic 
and fully ferromagnetic 2D HEGs, extrapolated to the ther- 
modynamic limit. 



Mag. state r 3 m* (a.u.) 



Para. 
Para. 
Para. 
Ferro. 
Ferro. 
Ferro. 



1 0.955(2) 
5 1.04(2) 

10 1.03(4) 
1 0.851(5) 
5 0.74(1) 

10 0.70(3) 



current data is much smaller, allowing a systematic trend 
with system size to be discerned and hence removed by 
extrapolation. 

For the paramagnetic HEG, GW calculation s 9 ' 28 indi- 
cate a steep increase in the effective mass as the density 
is lowered, similar to that seen in early experiments^ 
However, the GW results depend strongly on the choice 
of effective interaction and whether or not the calcula- 
tions are performed self-consistently. The QMC calcula- 
tions of Holzmann et al£2- give quite different results from 
those of either Kwon et al.^ our previous work^S or the 
present work. The experimental datai^ show some evi- 
dence for enhancement of the effective mass at low den- 
sity, although we do not see this in our present results. 

For the ferromagnetic case, our effective-mass data are 
in agreement with the experimental results of Padman- 
abhan et al.^ showing a decrease in the effective mass 
as the density is lowered. GW theory^ also predicts a 
suppression of the effective mass in the range of den- 
sities considered. However, the difference between the 
GW results obtained self-consistently and in the on-shell 
approximation is significant, as is the difference with the 
present results. 
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FIG. 6: (Color online) Quasiparticle effective mass m* against 
density parameter r 3 for paramagnetic or partially spin- 
polarized 2D HEGs, as calculated or measured by different 
authors. The experimental results are due to Smith and 
Stiles^, Tan et a/.,— and Padmanabhan et ai.— The GW re- 
sults were obtained using the random-phase-approximation 
(RPA) effective interaction^ and the Kukkonen-Overhauser 
(KO) effective interaction^ by solving the Dyson equation 
self-consistently (SC) or within the on-shell approximation 
(OSA). We show the VMC results of Kwon et a/.— (which 
were later confirmed at the same system size at r s = 1 a.u. us- 
ing transient-estimate DMC calculations^), the VMC results 
of Holzmann et aZ.,— and the DMC results reported in our 
previous work,— as well as the results of the present work. 
All the results shown are for paramagnetic HEGs with the 
exception of the experimental results of Ref. [J, which are for 
a partially spin-polarized HEG. 



CONCLUSIONS 



We have used DMC to calculate the single-particle en- 
ergy band and hence quasiparticle effective mass of the 
2D HEG. We have achieved sufficiently high precision in 
our calculations that systematic finite-size errors in the 
quasiparticle effective mass can be observed and removed 
by extrapolation. This leads to a revision of the effective 
masses for paramagnetic HEGs at low density compared 
to our earlier work)^ in particular we find that there is 
no enhancement of the effective mass at low density. 
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FIG. 7: (Color online) Quasiparticle effective mass m* against 
density parameter r 3 for ferromagnetic 2D HEGs. The 
GW results were obtained using the Kukkonen-Overhauser 
(KO) effective interaction by solving the Dyson equation 
self-consistently (SC) or within the on-shell approximation 
(OS A).— The experimental results are due to Padmanab- 
han et a/.— We show the DMC results reported in our earlier 
work^S in addition to our current results. 
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